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Abstract 

^ . We study the two-dimensional repulsive Hubbard model by functional RG methods, using our 

recently proposed channel decomposition of the interaction vertex. The main technical advance of 
this work is that we calculate the full Matsubara frequency dependence of the self-energy and the 
interaction vertex in the whole frequency range without simplifying assumptions on its functional 
form, and that the effects of the self-energy are fully taken into account in the equations for the flow 
of the two-body vertex function. At Van Hove filling, we find that the Fermi surface deformations 
remain small at fixed particle density and have a minor impact on the structure of the interaction 
vertex. The frequency dependence of the self-energy, however, turns out to be important, especially 
at a transition from ferromagnetism to d-wave superconductivity. We determine non-Fermi-liquid 
exponents at this transition point. 
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I. INTRODUCTION 



The phase diagrams of the cuprates and other currently investigated materials exhibit a 
variety of distinct phases, some with symmetry breaking, others with anomalous normal- 
state properties. Different effective theories for the competition of ordering tendencies have 
been proposed and compete with each other. On the atomic length scale, the band structure 
of these compounds has been determined and a description can be given by a three-band cor- 
related fermion model.-*^ By further arguments, one can reduce it to a single-band Hubbard 
model on a square lattice, with a Fermi surface consistent with ARPES measurements.-!^ 
In these models, the electronic Coulomb interaction is replaced by a short-range repulsion. 
In this paper we consider an on-site interaction. This class of itinerant fermion models was 
introduced a long time ago,^"^ but has received much attention in the context of high- Tc 
theory, and is of great theoretical interest by itself, due to the subtleties of the interplay of 
the interaction and hopping terms. Major aims are a controlled derivation of an effective 
low-energy model for the cuprates, and in the end, the calculation of the full correlation and 
vertex functions of these models. In particular, a full knowledge of the fermionic self-energy 
and the two-particle interaction vertex would allow to predict equilibrium and transport 
properties, to compare to Fermi liquid behavior in the normal phase, and to determine the 
dynamics of symmetry-broken phases. 

In this paper, we continue and conclude our analysis^'^ of the functional renormalization 
group (RG) flow in the level-2-truncation for the two-dimensional Hubbard model in the 
symmetric phase (all details will be specified below). We calculate the frequency and mo- 
mentum dependence of the fermionic self-energy and the two-body vertex function, taking 
the feedback of the self-energy on the vertex flow fully into account. For that purpose we 
solve the functional RG equations numerically on an adaptive frequency grid, hence obtain 
the frequency dependence without any simplifying assumptions on the functional form that 
it takes. We can therefore test the quality of ansatzes made previously, which in general 
turn out to be accurate only in a very limited frequency range, and we derive improved, 
yet still simple functional forms that are accurate for all but very large frequencies. We 
use this to map out the ordering tendencies in the Hubbard model and determine quasi- 
particle properties in the symmetric phase. At Van Hove filling, we thereby confirm earlier 
findingaiSi"— of a transition from ci-wave superconductivity to ferromagnetism as the ratio 
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of hopping amplitudes between nearest and next-to-nearest neighbors is varied. Our result 
is again consistent with a quantum critical point at this transition, and we determine the 
exponent for the frequency- dependence of the fermionic self-energy at this point. This ex- 
ponent is close, but not identical, to the exponent found in Ref. |l3|; the deviation may be 
due to the presence of the Van Hove singularity (VHS) on the FS in our case. Because the 
effective interaction remains small down to very low scales near the quantum critical point, 
our results also allow us to give reliable values for parameters of a detailed effective model 
in its vicinity. An analysis of this effective model using RG flows with bosons and fermions 
is in progress. 

We now briefly outline the background of the method and the context of our results. The 
functional RG gives a functional differential equation for the correlation functions as a func- 
tion of an energy or length scale, the RG scale. It is a controlled method at weak coupling — 
and has been instrumental in analyzing the competition of ordering tendencies,— 'i^'^^"— 
transport properties of systems with impurities^ as well as non-equilibrium phenomena.— 



For a recent review, see Ref. 
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In fermionic theories, the functional RG equation is equiv- 
alent to an infinite hierarchy of differential equations for the vertex functions. The main 
challenges in solving this hierarchy are (1) that the vertex functions have a nontrivial mo- 
mentum and frequency dependence, which is crucial for understanding the above-mentioned 
physical phenomena, and (2) that all ordering tendencies that turn up via continuous phase 
transitions are at some point associated with interactions that get long-range, hence the 
momentum- and frequency- dependent vertex functions get large and eventually singular at 
certain loci in momentum space as the RG scale is lowered. Our present work is a certain 
completion of addressing the first challenge in the scale regime where no long-range order 
has developed, namely we solve the RG flow equations without any further assumption on 
the frequency dependence, and with an already well-tested approximation for the momen- 
tum dependence. Obtaining an accurate frequency-momentum dependence is particularly 
important for quantitative calculations in the IPI scheme^^ where full propagators are al- 
ways supported at all scales above the RG scale; in fact, the f2-scheme^ employed here uses 
only a very mild regulator, for which the propagators are nonzero everywhere in momentum 
space and also at all nonzero frequencies. Schemes with too strong regulators, such as strict 
cutoffs, fail to capture ferromagnetic correlations.— The second major challenge mentioned 
above has been tackled in various papers using simpler approximations. To get accurate 



statements about the given microscopic model, both methods have to be combined, and we 
now discuss our approach to this. 

From the point of view of RG analysis, one can distinguish different regimes of energy 
scales - the "high-energy" regime of degrees of freedom far away from the Fermi surface 
(FS); the intermediate energies, where renormalization effects become sharper and order 
parameters start to emerge, i.e. selected couplings start to grow in the flow and self-energy 
effects become important; and the low-energy regime where symmetry breaking happens, 
gaps may open on (part of) the FS, and where order parameter averages and fluctuations 
play a central role. The distinction between these regimes can always be made conceptu- 
ally in weakly coupled systems, but it is particularly pronounced in very weakly coupled 
systems, where the non-perturbative symmetry-breaking phenomena generically occur at 
scales that are exponentially small in the inverse interaction strength, but self-energy effects 
arise already in low-order perturbation theory and may have important physical effects at 
temperatures above symmetry-breaking scales. In special parameter regimes, the growth of 
the flowing couplings may be suppressed altogether and hence the second regime can reach 
to very low scales; at quantum critical points, it can reach to scale zero. The fermionic 
RG has been used to access the symmetry-broken phase,— but the full dynamics of the 
Goldstone modes has been captured better in a representation with bosonic fields.— On 
the other hand, the fermionic RG has been more accurate at the higher and intermediate 
scales, and the ansatzes made in bosonic studies have indeed relied strongly on the results 
of previously done fermionic flows. 

In Ref. y, we have proposed a parametrization of the fermionic RG by a natural channel 
decomposition, which allows to keep most details of the frequency-momentum dependence 
and to switch to a bosonic description at a certain scale, using the semi-group property 
of the RG to stop the flow and a Hubbard- Stratonovich (HS) transformation to bosonize 
selected fermionic bilinears, namely those corresponding to the most important terms in 
the effective interaction. An advantage of this procedure is that the parametrization leaves 
little or no ambiguity in the HS transformation of the effective interaction. Moreover, the 
parametrization is numerically efficient, so that a large number of interaction terms can be 
kept, to test the accuracy. To avoid any bias in the flow, the bare on-site interaction was kept 
as a separate term in the parametrization of Ref. |8|. This was then also adopted in a study 
using boson fields.— The fermionic formulation furthermore has the advantage of allowing 



for terms that do not have a straightforward rewriting in terms of HS fields, hence may be 
missed in ansatzes for bosonic actions. One of the results reported in Ref. ^ are interaction 
terms of exactly this type, which were discovered only in the fermionic scheme and which 
cannot be neglected in the flow. More generally, a restriction to (anti)ferromagnetic and 
superconducting correlations does not describe the vertex function of the Hubbard model 
quantitatively in the filling regime interesting for the cuprates. 

Several RG studies including self-energy effects have been performed in the past. In the 
situation of a moving FS, one faces the problem that momentum space cutoffs around the 
free FS cannot provide an effective regularization. By putting a counterterm function in the 
quadratic part of the action and solving the emerging inversion problem, a map between the 
non-interacting and interacting FS was constructed and extensively studied^® for systems 
with a regular free FS. Alternatively, the scale decomposition of the propagator can be 
adjusted dynamically according to the moving FS.^^ With this method, the stationary self- 
energy was calculated for several momenta close to the (moving) FS and at fixed particle 
density in a parameter region with dominant pairing instability.^'' Third, the above FS 
cutoff problem can be circumvented by employing an alternative regularization. With the 
temperature flow scheme^ and in the ferromagnetic parameter region, the Fermi surface 
deformation was examined at zero and non-zero magnetic field by fitting the scale derivative 
of the frequency-independent self-energy for several momenta near the FS to an ansatz 
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show that the FS 



with two hopping correction terms.— The studies in Refs. 
deformation remains small and has a minor impact on the structure of the interaction 
vertex at the considered points in parameter space; the FS tends to become more fiat 
in the superconducting parameter region and to become more curved in the ferromagnetic 
parameter region. 

The self-energy flow equation in the standard IPI hierarchy combined with a stationary 
parametrization of the interaction vertex leads to the generation of a frequency-independent 
self-energy. In previous studies, an approximation to the frequency-dependence of the self- 
energy was obtained by calculating the projection of the vertex flow to zero frequency, but 
then inserting the integrated vertex flow equation (with the one-loop frequency dependence) 
into the self-energy flow equation. In this way, the self-energy gets a frequency-dependence 
as in two-loop "sunset diagrams" , in which the momentum dependence is that of the flowing 
vertex function. This approximation was used to calculate the flow of quasi-particle scat- 



tering rates,'^ quasi-particle weights,— — and to continue to real frequencies^^i^ (Ref. |32 



uses the Polchinski scheme, Ref. 
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the Wick ordered scheme). Parts of the self-energy 
can be given back to the flow of the vertex by assuming that the self-energy correction to 
the one-particle dispersion is of the same order of magnitude as the Z factor, and then 
incorporating this Z factor in the propagator.-^^-^-SA All these studies show an anisotropy in 
the quasi-particle properties with considerable effects close to the Van Hove points. In a 
setting with partial bosonization, the imaginary part of the frequency- dependent fermionic 
self-energy was accounted for through its values at the two lowest Matsubara frequencies 
and at the Van Hove point.— 

In the present work we use the RG scheme for irreducible vertices, employing the level-2- 
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as well as the vertex parametriza- 



truncation and Katanin replacement as described in Ref. 

tion of Ref. ^ with its extension to frequency- dependent vertices given in Ref. 19|. In a first 
calculation, the stationary self-energy is investigated by determining the flow of corrections 
to the hopping amplitudes and chemical potential of the free system. This allows for reso- 
lution of the moving FS and a study of its influence on the RG flow. We calculate the flow 
of the stationary self-energy at fixed particle density, chosen to be interacting VHP, in a 
wide range of parameter values. In a second calculation we study the frequency-dependent 
self-energy by discretization in frequency space. Because we use the channel decomposition 
of the interaction vertex we can go beyond the approximation where a frequency-dependent 
self-energy is constructed from a stationary vertex, and we determine the full frequency 
dependence of the imaginary part of the self-energy. 

II. FERMIONIC RG SETUP FOR THE HUBBARD MODEL 

We consider electrons on the two-dimensional discrete torus P = Z|, i.e. on a square 
lattice of sidelength L G N with periodic boundary conditions, and subject to the Hubbard 
Hamiltonian 

^ = J] ^(P) C+,Cp,, + J] cJ+k, + 4-k,-Cq,-Cp,+ . (1) 

per* ' ' k,p,qer* 

s 

Creation and annihilation operators Cp^i are associated with particles of momentum p e 
P* = (^Zl)^ and spin projection s G {±}. The tight-binding dispersion 

e{x^ y) = —2ti (cos x + cos y) + 4^2 (cos x cos y + 1) — (2) 



describes particles hopping between nearest and next-to-nearest lattice neighbors. We con- 
sider ratios < ti < i. Here a chemical potential is already included, and /i = 
corresponds to free VHF. The screened Coulomb interaction is mimicked by an on-site re- 
pulsion U > 0. The model exhibits a non-trivial interplay of the kinetic part, which is 
diagonal in momentum space, and the potential part, which is diagonal in position space. 

In the path integral formulation, the action 

Sq{iP,iP) = J dp^qn{p) ippsipps -U j dk dp dq ipp+k,+'ipq~k-i'q-i'p,+ ■ (3) 

is a function of Grassmann variables ipps, ipps that are labeled by frequency- momentum 

tuples p = {po, p) G Mn X r* and spin projection s. Fermionic Matsubara frequencies read 

Po = nT (1 — e~*'^^"/"')/'i with po ^ = {1, 3, . . . , 2n — 1} and result from the division of the 

interval [0, 1/T] in 72 G N time slices. We adopt shorthand notation J" &p — |^ X^po^A/n per* 
in 

and 5{p — q) = y-5po,'jo'^p,q- Subsequently, the zero temperature limit T — )■ will be taken. 

To do the RG, we use the f2-scheme, where the inverse free propagator is multiplied by 
a soft frequency regulator,- 

qniv) = (ipo - £p) / Xn(po), Xfi(w) = (4) 

The scale parameter 1] > sets a scale in that for ^ Q, the regularized propagator is 
very close to the original one, while the pole of qo{p)^^ at po = and p on the FS gets 
replaced by a value of qn{p)~^ that is of order Q. In the limit Q 00 the free propagator 
vanishes whereas for = the regularization is removed. This regulator, as opposed to a 
momentum shell cutoff, allows for a simple regularization with moving Fermi surface due 
to self-energy effects. Moreover, small-momentum particle-hole processes, which become 
important for t2 > 0.3ti and VHF, are not artificially suppressed.— 

We use the regularization parameter Q as the RG scale parameter, and calculate the 
fermionic two- and four-point functions within the RG scheme for one-particle irreducible 
(IPI) Green functions, employing the level-two-truncation.— We choose an ansatz for the 
IPI generating functional that preserves all symmetries (translational, spin SU(2), charge 
U(l)) of the action and restrict subsequent calculations to the symmetric regime. 



The general form of this symmetric IPI functional ia^ 

s 

- ^ I dpi...dp4 "Y S{pi + P2 - P3 - Pi) ^pi,Si^p2,S2^P3,^3^P4,S4 
S1...S4 

^(^n(Pl,P2,P3,P4)'5siS4'5.2S3 - Vn{Pl,P2,Pi,P'i)5srsz5s2s}j- (5) 

Here is the (thermal) self-energy, and vq, the two-body interaction vertex. We write 
■^^^(^15^25^35^4) as a function of four frequency-momentum arguments, but by translational 
symmetry it only depends on three of them. Subsequent equations should therefore be read 
in the subspace Pi + P2 — Pz — Pa = 0- 

The RG equations for the above coefficient functions read^^ 

^n{p) = ]^ j d/ sn{l) (vq{p,1,p,1) -2vq{p,1,1,p)^, 
vn{pi, ...,pi) = {Tpp + Tph,cr + Tph,d){pi, ■ ■ .5^4)5 (6) 

where 

TppiPi, ■ ■ ■ 5P4) = ~\ j Ln{l,pi +P2-I) vn{pi,p2,l,Pi +P2-I) vniPi+P2 - IJ,P3,P4), 
Tph,cr{pi, ■ ■ ■,Pi) = ^\ j Ln{l,Pi -P3 + l) vn{pi,l,P3,Pi -Pa + ^nipi - Ps + l,P2,l,Pi), 
TphAPu-'-^Pi) = +\ j Ln{l,p2-P3 + 1) (2vn{pi,P2-P3 + iJ,P'i) vn{l,P2,P3,P2 - P3 + 

- vniPuP2 -P3 + hUPi) vnil,P2,P2- P3 + l,P3) 

- vn{Pi,P2 -P3 + l,P'iJ) vn{l,p2,P3,P2-P3 + 0) (7) 

denote the particle-particle, crossed and direct particle-hole contribution, respectively. The 
so-called single-scale propagator is given by sn{p) = —gn{p)4n{p) and the full propagator is 

9n{p) = {qn{p) + ^n{p)y'. (8) 

We replace^ L^{pi,p2) = Sn{pi)gn{p2) + gniPi)sniP2) ^ ^ (gniPi) gn{P2)) ■ By this modi- 
fication of the IPl fiow equations the feed-back of parts of the irreducible six-point vertex to 
f Q is included. This reorganization of fiow equations was found to be crucial for recovering 
the exact solution to the reduced BCS model within the RG method.— 
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In the limit f2 — t- oo, the microscopic action ([3]) should be recovered. This is why 
we supplement the system of differential equations with initial conditions via imposing 
Vfi{pi, . . . ,^4) — !■ 2U and T,fi{p) — )• —U/2 as ^2 — )■ 00. The non-zero limit of the self-energy 
for vanishing propagator is a special property of the Q regularization, J dl qQ^{l) -^n~^oo 1/2. 

Once the flow equations have been derived it is convenient to take the limit of continuous 
time^ n — 7- 00. We work directly in the infinite- volume limit L — t- 00, and also restrict to 
zero temperature, hence take T — )■ at > 0. Frequency-momentum tuples then take values 
P = (Po? p) € K X [— TT, tt]^, the inverse free propagator becomes qn{p) = {ipo — £p) /xn{Po), 
and we have 6{p — q) = {2n)6{po — go) (27r)^(5(p — q) and J dp = J dpo J dp, with J dpo = 
(27r)-i 4 dpo and /dp = i2n)-' dp. 

At zero temperature and VHF, a flow to strong coupling is observed: The interaction 
vertex develops a strong frequency-momentum dependence with divergences for several 
frequency- momentum pairs at a non- vanishing scale i7 > 0. We stop integration of the 
flow when the interaction vertex exceeds a fixed multiple of the free bandwidth, i.e. at 
max Iff^l = 40ti. This defines a "stopping scale" fl^, which has the interpretation of an en- 
ergy scale where correlations of particle-particle or particle-hole pairs become important^. 
To continue the flow to lower scales, a non-symmetric parametrization of the IPI functional 
is needed. The tendency towards a specific ordering is, however, already visible from the 
growth of the symmetric interaction vertex. 



III. EXCHANGE PARAMETRIZATION FOR THE INTERACTION VERTEX 

The parametrization for the interaction vertex proposed in Ref. [s] is designed to cap- 
ture the most singular vertex structure in a systematic way. Although used there only 
for frequency-independent vertices, it allows to include both the singular momentum and 
frequency dependence of the interaction vertex. Here we account also for frequency depen- 
dences and use this parametrization to study the RG flow including self-energy effects. 

The exchange parametrization is set up as follows: The solution to the flow equation ([6]) 
for the interaction vertex can be written Vfi{pi ■ ■ -Pi) = 2U + (f^ + fp/^ cr + ^ph,d) (Pi • • -^4)) 

with = Tpp, Vph^cr ~ %h,cr, i'^h,d ~ 'Tph,d- 

The differential equation for has the following structure (equations for and v^^^^^ 
are treated in the same way): The integrand on the right-and side contains the product 
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of two propagators, which exhibit singularities for certain values of the loop frequency- 
momentum, i.e. / or P1+P2—I have zero frequency and zero energy. Depending on the external 
frequency- momenta pi, . . . ,p4, these two poles can coincide. For such a configuration the 
flow is strongly driven, and it can be expected that this generates the most singular vertex 
structure. 

The feed-in of the interaction vertex itself to the flow is essential because it multiplies this 
propagator product, in particular its values near the propagator poles are of importance. If 
one assumes that the feed-in of the interaction vertex does not dominate the above mecha- 
nism, a simple form for the most singular structure of the interaction vertex can be written 
down. 

Since the external variables pi, . . . ,^4 enter the propagator product only via the linear 
combination pi + p2 the most singular dependence of on its arguments can be expected 
through the transfer frequency-momentum pi + P2- The remaining dependence on exter- 
nal variables should not be neglected since it can produce an important modulation of this 
singular structure. However, a specific form for this remaining frequency-momentum depen- 
dence can be assumed, for which previous studies using a Fermi surface patching^'^^-^^'' or full 
Brillouin zone discretization^^ can serve as a guideline, and which can be tested subsequently. 

Interpretation of f^, fp^^r; "^pl.d ^e found by regarding the spin structure of the 
interaction vertex. It turns out that these functions are directly connected to different 
interaction channels of two fermions, namely to interacting Cooper pairs v^q := f^, spin 
interaction v^j := v^f^^^j. and charge interaction := 2Vp,^^^ — T34fp^ T34 is the exchange 
operator defined as T34t;(pi,p2,P3,P4) = f (pi,P2,P4,P3)- 

We continue by using the functions v^f, Vg^, v^, then 



Following the above reasoning, each channel is Fourier decomposed in the non-transfer 
momenta in a way such that basic vertex symmetries (see equations (l20l) below) are satisfied. 
The pairing channel e.g. is written 



vn{p 



...,P4) =2f/ + {Vsc + Vm + ^Vk){pi, . . . ,P4) + -Vm{Pi,P2,P4,P3)- (9) 



%c(Pb---)P4) =D {JP1+P2 



Pi +P2 

2 



-Pi, 



Pi +P2 

2 

Pi + P2 

2 



Ps) 



Pi + P2 

2 



Pi)/«( 




-Pi. 




2 



m,n 



(10) 
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The full frequency structure of the interaction vertex including frequency-dependent 
boson-fermion vertex functions was studied in Ref. [ol. Non-transfer frequency dependences 
were found to be of minor importance to the parametrization of the most singular struc- 
ture of the symmetric interaction vertex. For the examination of self-energy effects we thus 
disregard non-transfer frequencies here. 

This results in rewriting the three interaction channels 



Pi + P2 ^ T^n ^ , X . /Pi + P2 



v^iPi ■■■P,) = -J2 U^^^ - Pi) DlAPi + P2) /, 



Pi -P3 



P3j 



Vm{Pi ■■■Pa) = +Y^ /m(pi ^^~^) ^^mn(Pl " Ps) fn{V2 + ^ ^ 

m,n 

v'^iPl ...P,) = -Y. UPi + ^1(P2 - Ps) /n(P2 - (11) 

m,n 

in terms of exchange propagators {M^„, -D^„, K^^} and form factors {fm}- In the following, 
scale parameter indices will be dropped. 

The parametrization (fTT]) does not uniquely determine exchange propagators since it 
disregards non-transfer frequencies. Equations f|TT]) can be solved for exchange propagators 
after specification of a choice for non-transfer frequencies. We consider these two frequencies 
as a function of the transfer frequency, thus projecting the frequency space to a line. All 
frequency projections that map a bosonic frequency to a fermionic one and respect the vertex 
symmetries are admitted. 

At temperature T = 0, non-transfer frequencies can e.g. be chosen as half of the transfer 
frequency. 



i,o_i.o_Ea ' 

"■1— '^3~ 2 



l,0_PO i,0__PO ' 
K.j_— 2 A2— 2 



Dmn{p) = - j dkidkg fm{^ " ki) /„(| - kg) Vsc{ki,p- h, h,P - h) 

Mmnip) = + y 'Jkldk2 /m(kl - |) ^(^2 + ^) ^^^(^1, ^2, h - p, ^2 + P) 

Kmn{p) = - dkidk2 fmO^i + ^) /„(k2 - ^) t;i^(A;i, ^2, k2 -p,ki+ p) (12) 

In the frequency-independent setup for systems at VHF, major contributions in the flow 
with exchange propagators come from diagonal propagator elements, namely exchange prop- 
agators associated with an s-wave form factor /i(x, y) = 1 in all channels as well as a (i-wave 
form factor f2{x,y) = cosx — cosy in the superconducting channel.- A subsequent study^ 
shows, by comparison to a flow with discretization of the interaction vertex itself as a func- 
tion of three fermion momenta, that this decomposition captures well the interaction vertex 
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'II' 2 -(ly 

fc2 fca 

Figure 1. Graphical representation of the parametrization ([9]) and (|lip of the interaction vertex 
with an s-wave form factor fi{x,y) = 1 in the three interaction channels and a d-wave form factor 
f2{x,y) = COST — cosy in the pairing channel. 

structure for several Fermi surface geometries. It also indicates that, in the parameter region 
of strong competition between FM and d-SC ordering tendencies, the d-wave form factor 
could be slightly modified, possibly in a scale-dependent way, in order to resolve the full sin- 
gular structure of the interaction vertex. This modification of /2 is currently investigated. 

We restrict the vertex flow by keeping only the functions Mu, Ku, Du, D22- The RG 
equations then read 



S(p)=-(-2f/ + irn(0)) / dls{l) 



+ s{l-p) 



1 



^ j dlsil+p) [Ku{l)+SMnil)] 

(13) 



and 



Dn.m{p) = +^ 1 d/ L (-(/ - |), / + I) (^-D^miP), ^^^Y^) (^0, 1, p), 

Mii(p) = ^ - f + f) (+Mn(p), ~^^+^^~^ ) (/o,l,p), 



m = 1,2, 



-2D + 3M + K 



(^0,1, P)- 



The feed-back of the interaction vertex to the flow is given by 



Fi(A,5)(Zo,l,p) =2U + A 

+ Jdu (i?n(/o, u) + i?22(/o, u) /2(-| + 1 - |)/: 
F2{A,B){l,,\,p) =Am 

+ y"du5n(/o,u) /2(u-l), 



(14) 



^-2(-^ + i + f; 



(15) 
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and M22 = K22 = is understood. 

The dependence of functions Fi on their momentum arguments can be written exphcitly 
by the help of trigonometric identities, 



F,{A,B){lo,l,p) = 2U + A 

+ {822)2 (lo) (cos(/. - |) cos(/. + |) + cos(/, - |) cos(/, + |)) 
+ {B22), (lo) (sin(/. - |) sin(/. + |) + sin(/, - |) sin(/, + |)) 

- (^22)4 (^0) (cOs(/^ - y ) COs(/y + y ) + C0S(/^ + y ) COs(/y - ^) j , 

F2iAB){l,,\,p) =(^A+ {Bu),{lo)) /2(1), (16) 
where (B). (/q) = / du B{Iq, u) 5'j(u) and 

gi{x,y) = 1, 

£/2(x,y) = - [cos ■^ + cos -j , 

.X X y 

g4{x,y) = cos -cos-, 

^5(a;,i/) = ^(cosa; + cos?/). (17) 

We remark that by the symmetries (i) Fi{A, B){lo, —l,p) = Fi{A, B){lo,l,p) and (ii) 
L(^±(/-|),/ + |) = L(^(^±(/o-2|),l-f),/ + |) of the integrand in ([H]), loop momentum 
integration can be restricted to half of the Brillouin zone (BZ). 

The essential feed-in of the rf-wave pairing channel to s-wave channels can be reduced to a 
simple form, which produces an approximation to the function Fi: The exchange propagator 
F>22{po, p), as a function of momentum, generically exhibits a pronounced maximum at p = 
and then decays rapidly. Thus, the momentum integrals (.022)3 (-^22)4 — (-^22)2 be 
small compared to {022)2- Moreover, the (£'22)3 term in Fi is multiplied with a function that 
vanishes at 1± |p = (0, vr), i.e. is suppressed near the Van Hove points, which dominate the 
flow at VHF. Dropping the (£'22)3 and (£'22)4 — {^22)2 terms in Fi yields the approximation 

F,{A,B){loA,p)^2U + A +(Sii)i(/o) + (^22)2 (^o) /2(1 - |)/2(1 + |). (18) 

13 



We have verified in numerous situations that this approximation produces results of high 
accuracy. In some setups this approximation can substantially lighten numerical efforts, it 
will be applied for the examination of the frequency-dependent self-energy, sec. |Vll 



IV. SYMMETRY CONSIDERATIONS AND DISCRETIZATION PROCEDURE 

The flow equations constitute a system of differential equations for functions. Parametriza- 
tion of frequency-momentum dependence for vertex functions will reduce the system to a 
number of coupled ODEs that can be studied numerically. This parametrization is guided 
by symmetry considerations. 

Let S{po,p) = (— PO)P) denote the operator of time reflection and T(po,p) = {po,Tp) 
an operator of spatial transformation. Let beT=<( ^'');(^ ^ j,(°^)>asetof 



1/ yo-iy \io 

spatial reflections. The solution to the flow equations ([6]) with the initial condition stated 
above exhibits the following symmetries 



and 



Tj{Sp) = (time reversal symmetry) 

E(Tp) = S(p) V T G T (spatial reflection symmetries) (19) 

v{p2,Pi,P4.,P3) = v{pi,p2,P3,P4), (anti-symmetrization property) 

v{P3,P4.,Pi,P2) = v{pi,p2,P3,P4), (partlcle-holc symmetry) 



v{Spi, . . . , Sp4) = v{pi, . . . ,^4), (time reversal symmetry) 

v{Tpi, . . . , Tp4) = v{pi, . . . ,^4) V T G T (spatial reflection symmetries) . (20) 

All these symmetries are a consequence of the specific ansatz ([5]) for the effective action. The 
anti-symmetrization property directly results from the expansion ([5]). The other symmetries 
are valid for the initial condition and in addition compatible with the flow equation, hence 
hold during the flow. The spatial reflection symmetries follow from the symmetries of the 
free dispersion relation e{x,y). 

As a consequence, exchange propagators exhibit symmetries 



Bmm{Sp) = Bmm{p), (time reversal symmetry) 

BmmiTp) = Bmmip) V T G T (spatial reflection symmetries) (21) 
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(0,0) 



Figure 2. Illustration of momentum discretization for exchange propagators in 12 segments about 
p = (0,0) and p = (vr,7r). Small circles mark the representative momentum associated with each 
segment. For detection of incommensurate AFM ordering tendencies placement of data points 
along the coordinate axes is advisable. 

for each of the propagator variables B = D, M, K. Furthermore, the combination of 
anti-symmetrization property and particle-hole symmetry implies D*{p) = D{p), M^{p) = 
M{—p), K^{p) = K{—p). Thus, diagonal elements in the magnetic and scattering channel 
are real, Mmmip) , Kmmip) e M. 

We continue with the description of our discretization procedure. The momentum de- 
pendence of exchange propagators is approximated by step functions, as in Ref. ^. This 
means we decompose the Brillouin zone into sectors centered about momenta p = (0, 0) 
or p = (tt, tt) and designed such that momentum regions with strong variation of exchange 
propagators can be resolved in detail. Exchange propagators B{p), p = {po,p), are thus 
written 

B{p)= Yl B'{po)xs{p), (22) 

Segments S 

with xs the characteristic function of segment S, and B^ describes frequency dependence 
inside segment S. By symmetry one can restrict the discretization procedure to one eighth 
of the BZ, e.g. to momenta {{x,y) : < y < x < vr}, see Fig. [2l The actual choice 
of segments can be individually adapted to the type of exchange propagators for optimal 
resolution. A quadratic spacing in the radial direction allows a detailed study of the D22 
propagator. Linear spacing is used for all other exchange propagators. 

The frequency dependence of exchange propagators inside each momentum segment can 
be taken into account in several ways. The simplest possibility is to assume constant func- 
tions B^[pq) = B^(j)q = 0). The restriction to zero frequency has been widely made in the 
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past and focuses on the frequency value of most singular behavior. 

The suitability of more elaborate parametrizations which account for the decay of ex- 
change propagators in the frequency variable has been examined in Ref. [ol. Motivated 
from the small/large frequency asymptotics of RPA resummations it is natural to em- 
ploy the specific functional form of a Lorentz-like distribution, i.e. a rational function 
uj (Op + ih-p w + Cp w^)"^ with momentum-dependent parameters ap,6p,Cp G M. These 
three parameters can then be determined from the flow equation e.g. from the values of the 
exchange propagators and their first and second order frequency derivatives at zero transfer 
frequency. However, it seems that only small-frequency behavior is captured well by this 
procedure, which is insufficient for precise calculations. 

Here we do not make an assumption on the functional dependence on the transfer fre- 
quency variable, but instead discretize it on a frequency grid with logarithmic spacing 
and include ujq = 0. This choice allows us to study in detail the small frequency region where, 
especially at low scales fl, exchange propagators can vary strongly. By symmetry it suffices 
to take into account only non-negative frequency values. 

The discretization of self-energy is specified in the respective sections. 

V. FERMI SURFACE FLOW 

In this section, we discuss the detailed setup and the results of a flow in which the 
frequency- independent part of the self-energy is fully taken into account. This allows us to 
determine the RG flow of the Fermi line. 

A. Flow equations for hopping correction parameters 

The frequency-independent part of the self-energy satisfles 

S(po = 0,p)=i(-2f/ + irn(0)) jdl s{l) 

+ s(/o, 1 + p) [Bu{1) + i?22(/) f^il + P)) , (23) 

where the interaction vertex enters via the linear combination of exchange propagators 
Bii = ^{2D—3M—K)ii and B22 = -D22, and s{l) = sn{l) denotes the single-scale propagator. 
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Here we use the approximation of a frequency-independent interaction vertex and set all 
frequency variables to 0. By symmetry, S(po = 0, p) G M. 

Note that S(0, p) satisfies the symmetries (I19p of the free dispersion relation. This is why 
we parametrize the frequency-independent self-energy as a sum of corrections to hopping 
terms 

S(0,p) = /i,(p), Q eM, (24) 

with functions hi chosen from an ONS w.r.t. (/, (?) = / dp /(p)5'(p)- The first elements of 
this system are 

ho{x,y) = 1, 

hi (x, y) = cos X + cos y, 

h2{x,y) = 2 cosx cosy, (25) 

h^i^x, y) = cos 2x + cos 2y, 

h4{x, y) = v^(cos 2x cos ?/ + cosxcos2?/). 

We set up a flow at constant particle density, which fixes the coefficient Cq; details are given 
below. The remaining coefficients {ci)i>i can be determined in different ways, which in 
general will give different results. The first is a determination of the Cj from a global average 
in the sense: From the ansatz fl2^ the flow of coefficients Cj, z > 1, is uniquely given by 
orthogonal projection of il(0, p). Orthogonal projection Cj = (^hi, il(0, ■)^ yields 

c^ = -\ f dl s{l) R,{\), {l > 1), 

/ (26) 
i?,(I) = jdph,{l-p) (i?n(0,p) + 522(0, p)/|(l-|)). 

Here, the dependence of Ri{l) on 1 can be extracted from the momentum integral again by 
trigonometric identities. We have furthermore Ri{Tl) = Ri{l) VT G T, i > 1, since the 
functions hi have this property. 

Another method to flx the Cj is by local information at special points in momentum 
space, such as a Taylor expansion of S(0, p) around the Van Hove points. Comparison with 
the proposed parametrization then results in a linear system that can be solved for the q. 
Whereas the above orthogonal projection method allows to determine the q uniquely, the 
local expansion of S(0, p) in general involves all hopping modes that appear in (12^ . There- 
fore, a flnite linear system using local information extracts an approximation to (flnitely 
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many) coefficients q. Focusing on the saddle point region p = (0,7r) is natural since the 
Van Hove singularity substantially drives RG flows. In principle, the Hessian (9^11) (0, (0, tt)) 
provides two equations for flxing c\ and C2. However, using only two expansion parameters, 
we observe a discrepancy to the above orthogonal projection method. In fact, this highly 
local procedure completely disregards self-energy flow for momenta away from the saddle 
point region, which, though, signiflcantly feed in via the constant particle density condition, 
and can lead to unstable behavior. 

This is why we supplement the linear system by two further equations such that focus 
on the saddle point region is lifted. We consider the (non-degenerate) system 

((-a^ + a2)E)(0,(0,7r)) = 2ci 

((9.^ + 9.^) S) (0,(0, tt)) = 4c2 

E(0,(0,0))-E(0,(0,7r)) = 2ci +4c2 

S(0, (tt, tt)) - S(0, (0, tt)) = -2ci +4c2 

Below, we compare the results obtained from the two methods. 



C4 



-8c, 



+2V2C4 
-2V2C4 



(27) 



B. Fixing the particle density 

In an RG setup with flowing self-energy the Fermi line can change shape and level during 
flow. Consequently, the particle number will in general become scale-dependent. We consider 
the flow at constant particle density 

:= y d/c g^{k) = n, (28) 

i.e. we adjust the self-energy zero mode cq such that nn — 0. From the propagator relation 

g = s — g^T,, the condition of constant particle density n = reads 

0^/dp (^^0 - Sp) x(Po) ^ ^r^ ^Hpom _ (29) 



(^ipo - £p + x(po)S(p)) (^ipo - £p + x(po)S(p)^ 

and can easily be solved for cq, since Cq enters only linearly. 

We choose the particle density such that the Fermi line of the system at the stopping 
scale with effective dispersion :— £p — En^(0,p) contains the Van Hove points. We 
loosely refer to this situation as "interacting Van Hove fiUing" , although full removal of the 
regularization at T = requires a non-symmetric vertex parametrization. 
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Figure 3. Flow of hopping corrections in the orthogonal projection scheme(thick lines) vs. the local 
expansion scheme(thin lines). The dashed line corresponds to considering only the linear system 
from the second order derivative of S at the saddle point, the solid curve results from extending 
this system by 2 further equations as described in the text. Both schemes agree well except in a 
small parameter region about t2/ti = 0.35, where the local expansion method misses the flow to 
interacting VHF. All energy scales are given in units of ti. 



The adjustment of the particle number n is accomplished via proper choice of the chemical 
potential fi, i.e. several RG flows with different n need to be calculated, from which the 
correct one is selected. We choose the value of the chemical potential fi such that the free 
system has filling n and do not modify fi during flow, as this would leave the framework 
of IPI flow equations. The self-energy degree of freedom corresponding to the chemical 
potential is the zero mode Cq. Its initial value is chosen such that the interacting system at 
initial scale Qo has filling n. This is the appropriate initial condition in the flow at fixed 
particle density. 
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C. Numerical implementation and results at Van Hove filling 



In the numerical implementation of the flow equations, frequency-momentum loop in- 
tegrals are evaluated as follows: By our choice of the regulator function (jl]), and because 
we take only Sf7(0,p) into account in this section, the frequency integral can be done an- 
alytically using contour techniques. To do this for the loop integrands Lj^, which contain 
[ipo — £:p + Xn{Po) S(0, p)]^^, the zeros of polynomials of third degree have to be calculated. 
The remaining momentum integral is then evaluated numerically. All calculations have 
been performed for interaction parameter U = 3ti. Exchange propagators are discretized in 
momentum space by 60 segments per eighth of the BZ, see Fig. [2J 

We implement the flow at constant particle density, chosen to be interacting VHF. This 
condition flxes the flow of the coefficient cq and thereby essentially determines the flowing 
effective saddle-point level £^^{0, vr) = £:(0, vr) — Sn(0, (0, tt)) which is of central importance: 
it determines the distance of the Fermi level from the Van Hove points and thus triggers a 
strong enhancement of the flow of the interaction vertex. 

During each RG step, scale derivatives are determined in the following order: First, (ci)i>i 
are calculated from eq. (12B]) : second, eq. (1291) provides Cq; flnally, the scale derivative of the 
interaction vertex is computed from eqns. (fT4|) . 

The hopping corrections calculated from the RG flow remain very small compared with 
the initial parameters (thick line in Fig. [3]). During the flow, ti and t2 get a correction of 
at most a few percent. Figure H] shows the value of correction parameters Ci, . . . , C4 at the 
stopping scale and the corresponding flne-tuned particle density for interacting VHF. The 
dominant hopping correction from the stationary self-energy is between nearest neighbors, 
as given by the coefficient ci. When altering the free FS geometry, ci changes sign at 
^2/^1 ~ 0.38. As a consequence we observe that the effective interacting FS, determined 
from ^p*^, is less curved than the non-interacting one for ^2/1^1 ^ 0.4, for larger ratios it is 
more curved; this is consistent with previous results.— 1^ The adaptation of the FS patching 
scheme to the moving FS in combination with a certain assumption on the radial momentum 
dependences of the interaction vertex similarly finds a considerable impact of the self-energy 
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zero-mode whereas the change to hopping parameters (t2/^i)eff is not too strong.— Ref. 
gives the FS deformation as a function of the filling n at ^2 = 0, calculated from second-order 
perturbation theory. 
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Figure 4. Left: Correction parameters ci to C4 at the stopping scale, displayed by the solid, 
dashed, dotted, dash-dotted line, respectively. Right: Particle density corresponding to interacting 
VHF in the orthogonal projection scheme with 2 (solid) and 4 (dashed) hopping correction terms. 
The plot shows the difference to free VHF. 



The fact that S(0,p) remains small during the flow results in a minor modification of 
the structure of the interaction vertex when including the FS flow. In Fig. |5] we compare 
the stopping scale in the setup disregarding self-energy (dashed line) to the calculation done 
here taking account of self-energy corrections through coefficients ci to C2 and ci to C4 (two 
solid lines, almost identical). Consideration of the FS flow only slightly alters the stopping 
scale. Moreover, parameter regions corresponding to different ordering tendencies practically 
do not change. This indicates that the coefficients C3, C4 have almost no influence on the 
interaction vertex flow, hence that a parametrization of S(0,p) by the first hopping terms 
suffices. 

It is an interesting test for the parametrization whether the local expansion method 
f l271) can produce a good approximation to the flow of hopping corrections. Concerning the 
coefficients Ci and C2, we find that this is indeed true for many parameter values. It is 
essential here not to focus on the saddle point region only when setting up flow equations. 
In this check we employ the particle density value for interacting VHF identified within 
the orthogonal projection scheme. The two last equations of the system (12 7p then provide, 
compared to consideration of (9pS)(0, (0,7r)) only, corrections that bring Ci and C2 close to 
the scale derivative obtained by orthogonal projection, compare the thin solid and dashed 
blue Une in Fig.|3l The same is expected for the parameters C3, C4 after appropriate extension 
of the linear system. 
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Figure 5. Influence of the moving Fermi line on the flow of the interaction vertex for systems 
at VHF. Tracing S(0,p) leads to a small change in the stopping scale compared to neglecting 
self-energy (dashed line). We have calculated self-energy flows with 2 as well as with 4 hopping 
correction terms; this no more produces a difference in the stopping scale, the two solid lines are 
practically indistinguishable. Symbols indicate most dominant ordering tendency: commensurate 
AFM (filled square), incommensurate AFM (open square), d-SC (circle), FM (diamond). 



Only in the region t2 ~ 0.35ti of competing pairing and ferromagnetic ordering tendencies, 
differing flows for the first two hopping corrections are found. This is mainly due to the fact 
that deviations in ci and C2 result here, via the constant particle density condition f l29|) . in 
a significantly deviating flow of Cq such that interacting VHF is reached in the orthogonal 
projection scheme whereas it is not reached in the local scheme. This difference has a strong 
influence to the flow of the interaction vertex and leads to qualitatively different behavior. 
Given the strong influence of the effective saddle-point level, the flow is sensitive to which 
information is used to fix the particle density via fl29l) : determination of the particle density 
from an only locally known S(0, p) can be unreliable, and the orthogonal projection method 
should be preferred. 

Away from this parameter region and at interacting VHF, a parametrization of S(0,p) 
by the first few hopping terms provides consistent results. We observe that the FS flow has 
almost no impact on the interaction vertex. 
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VI. FREQUENCY-DEPENDENT SELF-ENERGY 



A. RG setup 

Within the stationary vertex approximation, the flow equation (jH]) generates a frequency- 
independent self-energy. A non-trivial flow equation for the frequency-dependent self- 
energy can be obtained though from the stationary interaction vertex by inserting its scale- 
integrated flow equation into the one of the two-point function.— i^i^i^ The parametrization 
of the frequency dependence of the interaction vertex enables us here to calculate the full 
frequency dependence of the self-energy, beyond the above approximation. To this end the 
set of RG equations for the frequency-dependent interaction vertex^ is extended by a self- 
energy flow equation. Furthermore, the full self-energy feed-back to the flow is accounted 
for by keeping full propagators on the right-hand side of flow equations. 

The most singular frequency dependence of self-energy is expected for its imaginary 
part and at small frequencies. Analysis of in second order perturbation theory for 
a Hubbard-like system at VHF and with a square Fermi surface shows^ a logarithmic 
singularity in the first frequency derivative for momentum p = (0,7r): At zero temperature 
and for frequencies < u < |, 

((9^ImE(2))(a;,p= (0, vr)) = Ci |ln w| V C2 |lna;| + (30) 



with Ci,C2 > 0. It is also shown in Ref. |43| that VS remains small to all orders in the 



coupling, even at VHF. Therefore, the growth fl30l) can have a drastic effect. 

Perturbation theory for the Hubbard system with Q regularization, at VHF and for a 
curved Fermi surface, < t2/ti < |, yields a similar picture. Numerical calculation of the 
Z factor, i.e. the inverse quasi-particle weight, 

Zp = l + (c'^lmS)(u; = 0,p) (31) 

to second order in U yields a logarithmic divergence (with power 1) in the scale parameter 
Q for p = (0, 7r) and indicates regular behavior for all other momenta, see Fig. [61 

For a detailed study of effects from the frequency- dependent self-energy we consider the 
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Figure 6. Second order perturbation theory for the Zp factor at VHF and momenta p = (0, vr), 
(xo,xo), (0,0), (vr, vr) (sohd, dash-dot, dashed, dotted hues, respectively; p = (0, 0), (vr, vr) almost 
coincide in the first two plots). The abscissa xq provides the intersection point of FS with BZ 
diagonal. 

flow equation for its imaginary part 



Im S {p) 



lms{l+p) {iReD + ^K + ^M)n(0 + ReD22(/)/|(^ + 

1 



+ Res{l+p) [ ImDuil) + Im 022(1) fi{- + 



(32) 



The propagator including the imaginary part of the self-energy reads g{p) = xi^) — Sp + 
ix{^) IniS(p)) ^ and enters the vertex flow equation (|T4l) via 



9{Pi)9{P2) 



X(wi) 



^i2U(^2)ImS(p2)) 



X{uji) x{^2) - £pi + «x(^i) Im S(pi) iu2 - e^^ + %x{^2) Im ^(pa) 

(33) 



Note that in the flow equations the full feed-back of the frequency- and momentum- 
dependent imaginary self-energy is kept. 

Since Im S is expected to show non-regular behavior only in the small frequency region, 
a linear frequency parametrization ImS(a;,p) ~ (Z^ — 1) a; appears natural. However, we 
find that this linear approximation leads to a strong artificial suppression of the interaction 
vertex during the flow, as compared to calculations using a discretization in frequency space. 
This strong suppression is due to the contributions from intermediate to large frequencies. 
Here the linear frequency approximation highly overestimates ImE(u;,p), which, in fact, 
decays to as |a;| — )■ 00. 
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B. Numerical implementation 



The flow equations ffT^ and (1521) . including ( 1551) and a similar equation for the single- 
scale propagator s{p), are studied numerically. For this purpose, exchange propagators are 
discretized in frequency and momentum arguments by 16 momentum segments per eighth 
of BZ, as described in eq. ( 122|) . Inside each momentum segment we use a logarithmic grid of 
length 20. .25 for resolution of transfer frequencies 10~^ < Po/ti < 10^. The smallest non-zero 
frequency value here is chosen depending on the expected stopping scale. Furthermore, the 
frequency value po = is included in the grid. The function u i— )■ lmS(u;,p) is discretized 
on a similar logarithmic frequency grid for a fixed set of momenta p. The most singular 
behavior of ImS is expected at the momentum p = (0, vr). Feed-back of ImS to the flow is 
done by spline interpolation of discrete data. All calculations are done for the interaction 
parameter value f/ = 3ti. 

Although the parametrization of ImS as a function of only one frequency-momentum 
argument is straightforward and easily accomplished in a discretization procedure with a 
comparably small number of parameters, the numerical evaluation of flow equations involving 
a frequency- and momentum-dependent self-energy turns out to be highly time-consuming. 
This is due to the presence of the full propagator G = (C^^ + S)^^ instead of the bare 
propagator C on the right-hand side. For the evaluation of the right-hand side integrals, we 
have set up an analytic loop frequency integration with contour techniques (i.e. after fit of 
the discrete frequency data to sums of Lorentzians) or with piecewise analytic integration 
of spline data^^ but found both to be impractical due to the high powers of the frequency 
variable in the resulting rational frequency integrand. Thus a numerical evaluation of the 
three nested integrals is necessary. 

In the parameter regions |t2Ai| < 0.2 and 0.45 < t2Ai < 0.5 we find stopping scales fl^ > 
10~^ti. Here, the full flow including feed-back of the frequency- and momentum-dependent 
imaginary self-energy has been numerically calculated. As the momentum dependence of 
Im p) turns out to vary not too strongly, a simple momentum set {(0, 0), (0, vr), (tt, tt)} 
is used for its discretisation. Results will be shown later in this section, see figs. [H [121 

For parameters 0.2 < ^2/^1 ^ 0.45, stopping scales fi* drop below lO^^ti. Because of 
numerical complexity we then proceed in a different way: For a fixed momentum value q 

25 



we decompose 

9[P) =7 X 7 ; 

lioj - ep + ix{^) ImS(w, q) j + ix{^) [ IniS(p) - ImE(cu, 

+ , J^'^'v , X(-). (34) 



a{p) a{p) {a + A){p) 
with a{p) = iu — Ep + ImS(u;, q) as well as A{p) = ^^(a;) ^ Im S(p) — Im 2(0;, q)^ . 



Now 



"(P) i (^w + Im q)) - £p 



(35) 



factorizes in a p independent term x{^) a term that depends on scale only through 
r{uj) = u + Xn{^) ImSn(a;, q). We insert this decomposition in the right-hand side of flow 
equations. For all terms that do not contain the function A, the loop integration can be 
reduced to the frequency integration by computing appropriate momentum integrals before 
the actual RG flow is started, similar to the treatment in Ref. [ol. 

To be speciflc, the flow equations for exchange propagators are of the form 

+ Xn{^-)Xn{^+) 



+ jd\ /[a, A](p,/,q) ^,,(p,I)^,,(p,l)| (36) 



with 



in - £i_p/2) {ir2 - ei+p/2) 
^ii(P,l) ^j2(P,l) 



in - ei-p/2y {ir2 - ei+p/2) 



~ ^0 and /[a, A] comprises all terms of the decomposition that depend also on 

A. 

Here ttj^(^o)V'i(P) 1) denotes the vertex feed-back to the flow as given by equations 
(fT6|) or (fT8|) . We choose to employ the approximation (fT8|) because it reduces numerical 
efforts; we have checked that it is of high accuracy in this setup by comparing speciflc scale 
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derivatives. Hence, functions ipjip,^) ^ {!> /2(1 — |p)/2(l + |p), /2(1)} are present in the 
j summations of eq. fl36|) . 

Momentum integrations F^.^^^^^H^.^^.^ can then be computed beforehand as they are 
independent of scale. This needs to be done for all discrete momenta {p} where exchange 
propagators are traced and for all combinations (V'iuV^j^) ^^at occur. The singularity of 
these momentum integrals along the "frequency" lines ri = or r2 = is essential to the 
flow and needs to be resolved in detail. The small-frequency asymptotics of the functions 
F, H can be calculated for the case ipj^ = ipj^ = 1 and for zero momentum exchange from 
the knowledge about the free Hubbard density of states. For example, at VHF and for 



with Ci, C2, C3 > 0. We use a logarithmic grid on both sides of the point ri = and similarly 
for r2, altogether a two-dimensional grid of 240 x 240 points. The minimal distance in the 
grid to the lines ri = or r2 = is carefully adjusted depending on the minimal value of 
the scale parameter Q that is needed during flow. The flow equation for Im S is treated in 
the same way. 

The loop integration for the terms including A remains a highly time-consuming compu- 
tational task because three nested integrations need to be performed numerically at every 
RG step and for each coupling. However, comparison in the parameter region with high 
stopping scales as well as further tests discussed below suggest that for the choice q = (0, vr) 
in the decomposition ( IMll . i.e. for the momentum value with expected strongest feed-in of 
Im S to the flow, these remaining terms containing A are well dominated by the first term. 

Thus, we can approximately calculate the flow by flrst neglecting these remaining terms. 
Then, as an error estimate, the approximated and full scale derivatives at the thus obtained 
stopping scale can be compared for speciflc couplings. Below, this approximation is discussed 
further. Notice that global momentum-independent feed-in of lmJ]{u), (0, vr)) still allows the 
determination of ImS(a;, p) at different points p in the BZ. 



< t2/ti < 



1 

2' 




(37) 
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C. Results at Van Hove filling 



Here we discuss results of the flow with frequency- dependent self-energy where the cal- 
culations have been performed using the approximation of momentum-independent feed-in 
of ImS(u;, (0,7r)) to the flow, as described above. For parameter values where the full flow 
could be calculated this approximation turns out to be very close to the true flow, details 
about this point are given in the next section. We compare the results to previous flows 
disregarding self-energy effects. 

The frequency-dependent self-energy feeds back to the flow via g(uj,p) = xi^) [ico{l + 
^ImS((X',p) — £p] ^. We show the function u ^ ^ImS(w,p) in Fig. [7] for an 

exemplary flow at different scales. The limit w — )■ of this function yields the momentum- 
dependent Zp factor. From this plot the region of validity of a linear frequency dependence 
of ImS can be seen: this region is roughly < Q. We have calculated ImE(a;,p) for 
p G {(0, 0), (0, tt), (tt, tt)} and find, as compared to results of second order perturbation 
theory, that Zp factors for all three momenta usually get enhanced in the RG calculation, 
cf. figs. Eland m In particular, 2'(o,7r) shows a strong divergence as ^2 — ?■ and hence becomes 
large at small scales. For t2/ti = 0.355, a value of 2'(o,7r) ~ 7 is reached at the stopping scale, 
whereas 2'(o,o) and ^'(^^7^) tend towards finite values in the flow. For parameter ratios ^2/^1 
close to but smaller than i, the 2'(o,o) factor is driven by an emerging singularity in the 
density of states, originating from the band minimum instead of the saddle point region. 
As a consequence, 2'(o,o) > ^(o.tt) during a substantial part of the flow, this effect is already 
present in perturbation theory. 

The fact that Zp factors become signiflcantly larger than one leads to a suppression of 
the flow of the interaction vertex, and especially at low scales the divergence of 2'(o,,r) has 
a strong impact. Thus, the frequency-dependent self-energy (over-) compensates the effect 
of enhanced divergences and rising stopping scales from taking the frequency dependence 
of the interaction vertex into account.- In Fig. [9] we compare stopping scales resulting from 
different vertex parametrizations. For small ratios ^2/^1 or ratios ^2/^1 near |, the flow of the 
frequency-dependent interaction vertex including Im S is quite close to a flow that neglects 
self-energy and frequency dependences at the same time. This has already been observed^ 
at )f:2 = 0, in a simpler approximation where only a Z factor is kept. For intermediate 
hopping ratios, the stopping scale is even further suppressed. In particular, in the parameter 
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Figure 7. Frequency dependence of Im T, for parameter values t2 = 0.355ti, U = 3ti. The function 
a; I— )• ^ Im T,{uj, (0, vr)) is shown for different values of the scale parameter. All energies are measured 
in units of ti. 




n/ti n/ti Q/h 



Figure 8. Flow of the Zp factors (inverse quasi-particle weights) for p = (0,7r), (0,0), (vr, vr) (solid, 
dashed, dotted lines, respectively) until reach of the varying stopping scales. The calculations for 
^2/^1 = 0.1 and 0.45 include feed-back of the full frequency and momentum dependence of ImS. 
For t2/ti = 0.355 the approximation of the momentum-independent feed-back of Im S((^, (0, vr)) is 
used. 



region of competing ferromagnetic and d-wave pairing ordering tendencies the stopping scale 
drops drastically to values beyond numerical resolution, for 0.34 < ^2/^1 ^ 0.35 we detect 
stopping scales ^^*/ti < 10~^. This behavior suggests a quantum critical point separating 
both ordering regimes, which was also proposed in previous studies.— 

Besides the lowering of the stopping scale we observe two additional differences to the 
setup with frequency- dependent vertex disregarding self-energy.- First, a region of dominant 
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d-wave pairing, which had disappeared in the calculation with frequency-dependent exchange 
propagators, is recovered. We see this change as a direct consequence of the small stopping 
scales that are obtained now, which allow effective generation of d-SC correlations driven 
from the magnetic channel. Furthermore, the scattering singularity, which had appeared 
in the study with frequency- dependent exchange propagators, is weakened and no more 
becomes dominant. However, it is still present. The right-hand plot in Fig. |9] shows the 
values of the most singular couplings in each channel at the stopping scale. The scattering 
coupling is at most of half the size of the most singular coupling. 

Thus, at VHF the RG flows (a) neglecting frequency dependences and self-energy effects 
on the one hand and (b) parametrizing frequency dependence of the interaction vertex and 
ImS on the other hand are qualitatively very similar. Certainly, they differ in the size of 
the stopping scale. However, both identify the same types of Fermi liquid instabilities and 
even detect them in almost identical parameter regions: for interaction parameter U = 3ti, 
we observe at hopping ratios 1^2/^1! ^ 0.2 a region of dominant AFM as a consequence of 
(approximate) nesting. Upon further increase of the FS curvature, relatively strong AFM 
correlations generate in a Kohn-Luttinger-like effect an instability in the pairing channel with 
d-wave symmetry. For ratios 0.35 < ^2/^1 < 0.5, formation of d-SC correlations abruptly 
stops due to emergence of a strong ferromagnetic instability, which crucially depends on the 
presence of Van Hove points on the FS. 

D. Non-Fermi-liquid frequency dependence of the self-energy 

In the parameter region ^2/^1 ~ 0.34 and for several ratios U/ti = 2.5 . . . 3.5, the competi- 
tion between pairing correlations with d-wave symmetry and ferromagnetic correlations de- 
celerates the flow of the interaction vertex. The regime of validity of the level-two-truncation 
is large and permits to trace the RG flow down to relatively small stopping scales. Previous 
works have suggested the existence of a quantum critical point in this parameter region.— 
The present findings, which take into account the frequency dependences of the interaction 
vertex and of the self-energy, support this scenario. 

The drop in the stopping scale goes along with a strong reduction of the quasi-particle 
weight at the Van Hove points during the flow. For momenta away from the saddle point 
region we detect a standard FL frequency dependence S(a;,p) — S(0,p) = {Zp — 1) iu + 
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Figure 9. Left: Comparison of stopping scales fi* at VHF for different vertex parametrizations. 
Top-down: (i) frequency-dependent vertex disregarding self-energy, (ii) frequency-independent ver- 
tex disregarding self-energy, (iii) frequency-dependent vertex including ImS(a;, (0,7r)). Dominant 
ordering tendencies: commensurate AFM (dashes, filled square), incommensurate AFM (dotted, 
open square), d-SC (solid, circle), scattering instability (dash-dot, triangle), FM (short dashes, di- 
amond). Right: Values of largest couplings in the different channels, for setup (iii), at the stopping 
scale. 

0{uj In" |a;|) as w — )■ 0, with Zp reaching a finite value as — )■ 0, see the mid plot in Fig. [HI 
On the other hand, Im E(co', (0, vr)) gets substantial contributions from the flow at all 
scales Q: it grows strongly in the frequency region < but changes little for larger 
\uj\ /Q. As a consequence, the curves u i— )■ ImS(co', (0,7r)) at different scales Q coincide for 
frequencies 1^1 > lOQ, see Fig. HDl hence provide access to the frequency dependence of the 
limit i7 — 7- 0. In particular, the RG calculation gives the frequency asymptotics 



Im S((X', (0, tt)) ~ const |u;|"sgnco'. 



(38) 



with exponent a ~ 0.74 for frequencies < lO^^ti. This value of the exponent is consis- 
tently extracted at all scales Q < lO^^ti. The contribution (138!) dominates the iu term in 
the propagator at small frequencies and leads to non-Fermi-liquid behavior. 

Fermion systems at criticality are extensively investigated in the literature. In the case of 
a regular FS with constant non-vanishing Fermi velocity and interaction via a single bosonic 
channel, the Eliashberg resummation produces the non-FL frequency asymptotics (l38ll with 
exponent a = 2/3. However, for fermionic spin SU(2) interactions the Eliashberg expansion 
is unstable with respect to vertex corrections and self-energy corrections.— In the present 
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Figure 10. Small-frequency asymptotics of the self-energy at the Van Hove 

points. The discrete data shows ImS at different stages of the flow, = 
10-^ 3•10-^ 10"^, 3-10-^, 10-^ 3-10-3, 10-2, 10-1 (top-down). In the frequency 
regime lOfi < \uj\ <^ ti all curves coincide. This allows us to extract a (non-Fermi-liquid) 
frequency asymptotics ImS(a;, (0, vr)) ~ const |cj|°sgncj with a = 0.74 (dashed line). For 
comparison, the dotted line shows the asymptotics with a = 2/3, which is not met by the present 
data. 




situation of fermions with an aspherical FS including Van Hove points and with competing 
interaction channels, the functional form (l38|l remains, but we obtain a different exponent. 

In the same way, the RG flow provides information about the asymptotic frequency 
behavior in the various interaction channels. The identification of asymptotic regimes based 
on data for scales Q/ti > IQ-^ is less obvious in the case of the four-point function due 
to non-monotonic dependences on the scale parameter Q in the Q scheme: The exchange 
function Mu^oj, 0) increases during the flow for frequencies \u\ < fi, it decreases for \u\ > 
see Fig. [TTJ Assuming that the asymptotic frequency regime is already reached, we extract 
Mii(cj,0) ~ |i^r"^ with an exponent 02 ~ 0.13. . .0.2 (green line in Fig. [TTi) . Concerning 
the (i-wave pairing channel, it is hard to identify an asymptotic frequency regime from the 
RG data. 
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Figure 11. RG flow of the ferromagnetic exchange function oo i— )• Mn (w, 0) at different stages of the 
flow, = 10"^ 3•10-^ 10"^, 10-3, 10^2, 10-1, 10° (top-down viewed from w/ti = IQ-^). The 
dashed line ~ I'^l""^ represents the estimate for the asymptotic frequency behavior of Mii(cj,0) 
with the exponent 02 in the region 0.13 .. . 0.2, the line uses 02 = 0.13. 

E. Comparison to flow with momentum- and frequency-dependent self-energy 

Because of numerical complexity the flow with momentum-dependent Im S was calculated 
only for parameter values < ^2/^1 ^ 0.15 and 0.45 < ^2/^1 < 0.5. Here, the approximation 
of feed-in of momentum-independent ImS(ci;, (0,7r)) to the flow turns out to produce only 
a small error: Figure WI\ compares the flow of the dominant coupling for t2/^i = 0.15 (AFM 
regime) and 0.45 (FM regime), only a small correction is found when accounting for the 
momentum dependence of ImS. 

For the remaining parameter range 0.15 < ^2/^1 < 0.45, where flows go down to small 
scales, a larger error is expected. Although the full flow could not be obtained in this region, 
full scale derivatives for speciflc couplings can still be calculated here. This serves as an error 
estimate to the approximated flow with momentum-independent but frequency- dependent 
feed-back of ImS(ci;, q), q = (0,7r). Using the data obtained in this way, we compare full 
and approximated scale derivatives for the most dominant couplings at the stopping scale. 
Since ((9^ImS)(a; = 0,p) typically grows strongest for p = (0,7r) the approximated flow is 
usually more suppressed than the true flow. Since |ImS(u;,p) — ImS(a;, (0, 7r))| typically 
increases during the flow, the comparison of scale derivatives at the end of the flow yields 
the largest error in the scale derivatives during the whole flow. 
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Figure 12. Flow of the dominant instability in the AFM (left) and FM regime (right) in differ- 
ent calculation schemes. We compare the flow including ImS(a;,p) (solid and dashed lines) to 
several setups disregarding self-energy effects: frequency-independent vertex functions (dash-dot) 
and frequency-dependent vertex functions (dots). The solid line corresponds to the full flow with 
feed-in of the momentum- and frequency-dependent Im$](a;,p), the dashed one to a feed-in of 
ImS(u;, (0, vr)) only. At t2/ti = 0.45, ^(o,o) > ^{o,tt) most of the flow, which leads to the atyp- 
ical situation that the momentum-independent feed-back of ImE(a;, (0,7r)) slightly overestimates 
the true flow, instead of underestimating it. 

For parameter values t2Ai = 0.15, 0.3, 0.33 and 0.42 we find that scale derivatives of the 
leading and next-to-leading couplings of the interaction vertex get, at the stopping scale and 
by feed-in of momentum-dependent ImS, a relative correction of 1%..15%. 

Alternative choices q = (0, 0) or q = (tt, tt) underestimate Im S in the saddle point region 
and lead to larger stopping scales, see Fig. [131 With such choices of q, scale derivatives are 
highly overestimated, by a factor of 2.. 5, which underlines the special role of the saddle point 
region for systems at VHF. 

From this check it can be deduced that (a) accounting for the correct Im S in the saddle 
point region is of major importance at VHF and (b) the true flow with full Im E should be 
close to the one with Imll(ci;, (0, vr)) feed-back. 



F. Effect of the imaginary part of exchange propagators 

Already in perturbation theory a singularity in the imaginary part of the interaction 
vertex is found: Although the imaginary part of the particle-particle bubble vanishes at 
zero frequency, it shows a singular frequency derivative. For zero momentum transfer, at 
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Figure 13. Approximation of momentum-independent feed-back of ImS((^,q) to the flow. The 
results for several choices of q (solid: q = (0, vr), dashed: q = (0,0), dots: q = (vr,7r)) share 
three common characteristics compared to neglect of self-energy (uppermost curve): The stopping 
scale gets (for intermediate t2/^i substantially) suppressed; a d-SC dominated parameter region is 
recovered; the scattering process with non-zero frequency exchange is weakened. 
The flow with full ImS feed- in is expected to be close to the approximate one with q = (0, vr), 
but with a higher stopping scale. On the other hand, from the discussion in the text the true 
stopping scale is supposed to be well below the ones obtained from choices q = (0,0) or (vr, vr). 
Color conventions as in Fig. [9j 



This singularity is induced by the asymmetry of the density of states N'{e) about energy e = 
0. Upon ladder resummation in the SC channel, an even stronger singularity ~ (^7 — ^7*)^^ 
builds up in this frequency derivative when the composite pair field becomes massless. 

In our parametrization, only Dn and D22 can acquire an imaginary part, see the sym- 
metry discussion below eq. f l2T]) . We have already examined the influence of this singularity 
on the flow of the frequency-dependent interaction vertex with neglect of self-energy-S at the 
data point ti = 0.3, U/ti = 3. This influence was found to be minimal: despite the sin- 
gular frequency derivative, imaginary exchange propagators remain small during flow. For 
the same data point, we study here the influence of ImDu, ImD22 on the combined flow 
of the frequency-dependent interaction vertex and ImS(a;, (0, tt)). In this parameter region 



VHF and with < t2Ai < 
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Figure 14. Imaginary part of the D22 exchange function. Left: Flow of its first frequency derivative. 
Right: The function w 1— )• ;^ ImL>22(w, 0) at the stopping scale 0* w 0.002ti. 




Figure 15. Comparison of flows taking into account ImDn and ImZ)22 or neglecting them. Left: 
Flow of the leading couplings, line conventions as in Fig. [9l Right: Flow of -Z^(o,7r)- Both flows are 
almost indistinguishable. 

of dominant d-SC, AFM correlations are found to generate a strong ImD22)('^ = 0,0), 
which becomes much greater than the frequency derivative of ImZ^n, see Fig. [TH Never- 
theless the functions Im Dmm{^) remain small during flow, as found before. 

Accordingly, these imaginary vertex contributions practically do not change the flow of 
Im S and the most singular couplings, see Fig. [151 We leave the calculation of the flow with 
Im Dmm in other parameter regions for future work but do not expect a strong effect there 
either. 
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VII. FUNCTIONAL FORM OF FREQUENCY DEPENDENCES IN THE SYM- 
METRIC PHASE 



The functional form of the frequency dependence of exchange propagators and of the self- 
energy is of particular interest in view of going beyond the numerically expensive frequency 
discretization used so far. Based on the discrete RG data we now propose and discuss several 
parametrizations. 

RPA calculations motivate a transfer frequency dependence in form of a Lorentz curve,- 
such that the small-frequency region is captured accurately. However, we find that in the 
intermediate to large frequency region |ci;| > the true frequency behavior deviates from 
simple Lorentzians, especially at low RG scales, and that the intermediate frequency region 
l^l ^ Q can be of considerable importance for the calculation of a quantitatively correct 
flow. 

Generalizations of single Lorentz curves considerably extend the regime of validity of the 
functional parametrization for transfer frequency dependences u t— t- Bm,n{^,p) and for the 
frequency dependence of the self-energy in the form u t— t- Im(ci;, p): At not too low RG 
scales Q, a sum of two Lorentzians 

^(-' - ttI^ + ttI^ ''"^ 

captures frequency dependences quite well.-=S- These two curves have the interpretation of a 
small- and a large-frequency process. Furthermore, in the frequency region l^l < lOQ the fit 
to two Lorentzians is very accurate at all scales fl/ti > 10~^; this frequency region is most 
important in the calculation of the RG flow of the leading and sub-leading couplings. 

Whereas the exchange functions Mn, Dn and D22 in the magnetic and pairing channel 
remain positive for all frequency-momenta such that the corresponding fermionic interac- 
tion can be decoupled by a Hubbard- Stratonovich transformation, the RG flow detects a 
sign change in the scattering channel around zero transfer momentum with a maximum of 
Kn{u,0) at u; = and a pronounced minimum i^ii(c<;s,0) < at Ws > 0, see Ref. y. In 
the parametrization with Lorentzians, at least two Lorentz curves are needed to account for 
this behavior. 

The parameter range ^2/^1 ~ 0.34 exhibits very low stopping scales, and the ferromagnetic 
exchange function as well as the self-energy at the Van Hove points develop a characteristic 
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Figure 16. Overview about the various functional parametrisations for the transfer frequency 
dependence of exchange propagators at low RG scales. On the example of the ferromagnetic 
exchange oj i— )■ Mii(a;,0) several ansatzes for the discrete RG data (+ symbols) are compared: A 
single Lorentzian determined from the small-frequency behavior (red dashed line) is appropriate 
for small frequencies <^ 17. Likewise, the sum of two Lorentzians (green line in dashes and 
dots) is a compromise to the different behavior in different frequency regimes and not an accurate 
parametrization for a large frequency range. The functional form (|4ip leaves the frequency exponent 
in the region of intermediate frequencies free and describes the RG data quite well (blue solid line). 
It decomposes into the two blue dotted lines. The parameter values are ^2/^1 = 0.341, U/ti = 3. 

frequency asymptotics at small frequencies, see Sec. IVI D[ In this case, different frequency 
regimes occur: In the small-frequency region l^l ^ Q, the flow in the Q scheme detects at 
all scales a quadratic dependence as e.g. in Lorentzians; this result possibly depends on the 
regularization method. On the other hand, a pronounced intermediate frequency region with 
a characteristic decay \uj\~^ forms during flow, consider e.g. the region 10^^ < u/ti < 10^^ 
in Fig. [161 In the Q scheme, this decay behavior no more changes at decreasing Q for 
frequencies l^l > lOQ, and we therefore expect it to be independent of the regularization 
and to provide the small-frequency asymptotics as the regularization is fully removed. 

In this context, a suitable modification of fl40p is an ansatz where the frequency exponent 
7 > describing the intermediate frequency region can be varied. 



Fig. [16] shows the relatively high quality that is achieved by this ansatz, as compared to the 




(41) 
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parametrizations with one or two Lorentzians. All fits are determined with a least-squares 
method and a ten times larger weighting factor for frequencies \u\ < Q. The frequency data 
obtained in the RG flow has < 7 < 1 at low RG scales Q. A factor Xn{^) is included in 
the second term of fHT]) in order to satisfy f'{u = 0) = 0; this condition is imposed on the 
exchange propagators and on the function ImS(u;, p) by the Q regularization at scales 

> 0. For simplicity, the Q regulator Xni^^) = u)'^/{uj'^ + ^^^) is used in the ansatz. Apart 
from the scattering exchange Kii{p), which has a sign change along the frequency axis, we 
find that exchange propagators at all scales fi/ti > 10^^ are well captured by fHTl) . 

The frequency-dependent self-energy in the form u t— t- u^^ ImS(a;, p) is described well at 
all scales > 10^^ by the generahzation of fHT]) to 

l + bfu^ 1 + 02 \uj\ uj^ + 
where Xc{^) replaces Xn(i^) and c > is a further variational parameter. The ansatz 

/(-) = jmhw^ («' 

is designed to capture the quadratic dependence at small frequencies as well as the charac- 
teristic decay at intermediate frequencies with just three parameters. Whereas this ansatz 
does not reproduce the various transfer frequency dependencies, we find that (143|) indeed 
describes Im S very well in the small to intermediate frequency region and hence is par- 
ticularly suitable to account for self-energy effects when calculating the flow of the leading 
exchange couplings. Figures [16] and [T7| show the quality of these ansatzes. Tables [T] and [ITl 
provide the corresponding fit parameters. 

Notice that the fit parameters 7 are related to the exponents identified in Sec. IVI Dl the 
parameter 7 in the ferromagnetic exchange propagator corresponds to the exponent 02 given 
previously; the self-energy parameter 7 corresponds to 1 — a, since the fit uses an additional 
factor u~^. The fit parameters 7 in Tab. [T] are very close to the exponents determined in 
Sec. IVIDl however do not always coincide exactly. This is because 7 here is determined as 
the best compromise to describe the behavior for all frequencies, whereas the true exponent 
characteristic to the asymptotic behavior at small frequencies must be extracted exclusively 
from the frequency region where the asymptotics shows. 

It is interesting to note that the flow of the leading couplings (including their vicinity in 
frequency-momentum space) and the flow of the self-energy at small frequencies is in many 
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Figure 17. Effective functional forms from the discretized frequency dependences obtained with 
the RG flow for to i— )• uj^^ Im T,{io, (0, vr)) and to i— )• Re D22{(^, 0). The line conventions and further 
details are given in Fig. H The fit for Im S uses the additional freedom provided by the ansatz 
(j42|) instead of (j41|) : the cyan dash-dot-dot line shows the fit to the three-parameter ansatz (j43|) 
and captures well the discrete data except for the large-frequency asymptotics. The corresponding 
fit parameters are given in Tables [I] and |lTl Data at frequency a; = (not shown in the logarithmic 
plots) is very close to the values at the lowest frequencies given in the graphs, in agreement with 
the Lorentzian form at very small frequencies. The parameter values are t2/ti = 0.341, U/ti = 3. 
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Table I. Fitting parameters for the discretized frequency dependence of the ferromagnetic exchange 
Mii(a;,0) and the d-SC exchange D22(w,0), as obtained in the RG flow. At the considered pa- 
rameter point t2/ti = 0.341, U /ti = 3 these exchange channels are the most important ones. The 
fits to the ansatz ()4ip are provided at four different RG scales = 10^2 , , , 10^5_ Figures [T6l 
and 1171 show the quality of the fits. All the data in the table uses units where ti = 1. 
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Table II. Fitting parameters for the discretized self-energy frequency dependence w i— 
o;-^ Im (0, vr)), as obtained in the RG flow. The table gives the parameters for the ansatz 
([I3|) at several RG scales, the quality of this fit is shown in Fig. [T71 All the data in the table uses 
units where ti = 1. 



cases sensitive only to proper parametrization of a small frequency region < Q, especially 
at not too low RG scales. The behavior of exchange propagators and the self-energy in this 
small frequency region can be described well with a sum of two Lorentzians or, at not too low 
RG scales, even by a single Lorentz curve. The more general ansatzes (1411) to fj43l) capture 
well the behavior in a much larger frequency region at all scales fi/ti > IQ-^ and allow to 
describe the forming of the specific asymptotic frequency dependence as — )■ 0. 

The generalization of the functional form of the frequency dependence from a single 
Lorentz curve to (real- valued) functions ( l40l) - ( H3l) raises the question of how to extract the 
free parameters in the ansatz from the flow equation: while in these ansatzes two parameters 
can in principle be fixed by the Taylor expansion at zero frequency, the remaining parameters 
account for the behavior in the intermediate-frequency region. In the present work they are 
determined from a least-squares fit of discrete data. A more general way of extracting a flow 
equation for them needs yet to be established. With an appropriate method for setting up 
the flow equations for these parameters, the numerical complexity of the evaluation of the 
flow equations ([6]) will be substantially lowered and further insight in the structure of their 
solution can be gained. 
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VIII. CONCLUSIONS 



In this work we have solved the level-2 truncated RG equations for the two-dimensional 
repulsive Hubbard model, including the feed-back of the self-energy to the vertex flow equa- 
tions via the full propagators in the flow equations. We have computed the stationary 
self-energy and thus determined the band function of the interacting system, and, in a sec- 
ond step, calculated the imaginary part of the self-energy. To our knowledge this is the 
first time in this two-dimensional model that the frequency dependence of vertex functions 
and self-energy have been taken into account without simplifying assumptions on their func- 
tional form, while also keeping the momentum dependence in a well-tested approximation. 
Although our method is applicable in general, i.e. without restriction to particular parame- 
ter regimes, we have focused on the situation where the interacting system is at Van Hove 
filling, because of the intrinsic theoretical interest of this situation and because of the indi- 
cations from ARPES measurements that the Fermi surfaces of cuprates are indeed close to 
Van Hove points. 

The stationary self-energy determines the band dispersion of the electrons in the inter- 
acting system. In the RG setup, it becomes a scale-dependent quantity, hence the Fermi 
surface also fiows under the RG. This is a power counting relevant effect at weak interac- 
tions, hence a priori important. Capturing this fiow is technically nontrivial in momentum 
cutoff schemes, but much easier in the Q-scheme we use here since this scheme does not 
involve a momentum space cutoff. We have Fourier expanded the self-energy into a sum of 
correction terms to free electron hopping and fixed the particle density during the flow by 
adjusting the self-energy zero mode. By proper choice of the chemical potential, we tuned 
the density to interacting Van Hove filling. 

The corrections to the hopping amplitudes were calculated by orthogonal projection of the 
self-energy to standard hopping functions. This method was found to be consistent with an 
alternative one tliat extracts cocflicicnts from a Taylor expansion of the self-energy around 
the Van Hove points, supplemented by conditions obtained from the self-energy at other 
points in momentum space. We also found that the attempt to use only local information 
at the Van Hove points leads to unstable behavior, i.e. including higher terms in the Taylor 
expansion around these points leads to big changes in the coefficients determined from the 
lower orders. This problem disappears when an overemphasis on the Van Hove points is 
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avoided by taking into account other momentum space points. 

For the parameter values considered here and with our fixed-density constraint, we found 
only a small change of the Fermi surface in the fiow, resulting only in minor modifications 
of the flow of the interaction vertex. 

Our results on the imaginary part of the self-energy allow us to study several aspects of 
the frequency-dependent two-point function, most importantly the decrease of quasi-particle 
weights in the interacting system. We find that away from frequencies below or around the 
RG scale O, the hnear frequency parametrization of the self-energy, as obtained by Taylor 
expansion around frequency zero, is inappropriate and leads to an artificial suppression of 
the flow. Wc therefore discretize the sclf-cncrgy behavior in frequency space, putting special 
emphasis on the small-frequency region. In momentum space, the neighborhood of the Van 
Hove points is particularly important: perturbation theory finds a singularity in the first 
frequency derivative of the self-energy, and the RG calculation gives an enhancement of that 
singularity. This momentum space region also substantially drives the fiow of the interaction 
vertex, and hence properly taking into account the self-energy here is crucial. 

The above-mentioned problems with the accuracy of Taylor expansions are not surprising, 
especially for self-energies with a singular frequency behavior: Taylor expansion always 
works at very small scales, i.e. when \uj/fl\ is very small; this is in the nature of a regularized 
theory, where singularities are smoothed out by the regulator. It is, however, another matter 
to extend this to larger frequencies > Q, and our results show that the above-mentioned 
Taylor expansion around uu — does not represent the function in the frequency range above 
fl correctly. 

We have used our numerical results to fit specific simple functional forms to the self-energy 
and vertex functions. Our ansatzes provide a considerably larger regime of validity than 
simple Lorentzians coming from Taylor expansion, and will be of interest in future studies of 
the model, as they allow to avoid a fully numerical study while retaining the accuracy of the 
frequency dependence. The coefficients can again be determined with much less numerical 
effort (but - again - not only by a local Taylor expansion that would overemphasize the 
vicinity of the singularity) . The behavior at very large frequencies is given by a convergent 
expansion in inverse powers of ou, and may have coefficients differing from those we get in 
the intermediate frequency regime. This leads, however, only to small changes. 

Concerning the apphcation to the Hubbard model, we find that at the Van Hove filling, 
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as a consequence of decreasing quasi-particle weights, the flow of the frequency-dependent 
interaction vertex gets slowed down. This effect is most drastic in the parameter region 
of competing rf-wave pairing and ferromagnetic instabihties, where the stopping scale of 
the flow drops by several orders of magnitude, conflrming earlier suggestions of a quantum 
critical point. At this point we calculate the non- Fermi-liquid frequency dependence of the 
symmetric self-energy. It is a fractional power law, as found previously in models without 
Van Hove singularities, but with a different exponent. 

In comparison to the results from flows with a frequency- dependent vertex function 
where the self-energy is neglected,- the stopping scale for the flow of the frequency- 
dependent interaction vertex with self-energy feedback drops below that obtained in the 
static approximation .-li^ Similarly, a region of dominant (i-wave pairing reappears also at 
U = 3ti, and scattering processes with non- zero frequency exchange get weakened, consis- 
tent with the fact that the d-wave correlations also get driven by contributions from the 
zone diagonal, while the magnetic ones depend much more strongly on the behavior close 
to the Van Hove points. 

In summary, as far as the stopping scale and dominant correlations are concerned, the 
fully frequency- dependent flow calculated here agrees well with the results from our flows 
with frequency-independent functions. This is a reassuring indication that the projection 
to zero frequency, which gives the leading behavior in weak-coupling power counting, really 
works at the values of U that we consider. Our detailed analysis reveals, however, that this 
is not to be understood via a trivial Taylor-expansion argument, but that a more subtle 
mechanism is at work, namely that the decrease of the quasi-particle weights works against 
the relative enhancement of the vertex at zero frequency, which in turn is caused by the 
decay of the vertex function in frequency. 

The deviations appearing in the approximation where the vertex frequency dependence 
is kept but the self-energy is dropped are not surprising in hindsight; vertex function and 
self-energy are linked by a Ward identity and a fleld equation. It is more plausible that 
one can approximately keep these relations either by dropping the frequency dependence 
of both functions, or by keeping it for both, than by going only half-way. At this level of 
generality, this argument does, however, not explain why the static approximation also flts 
quantitatively for the above-mentioned quantities. 

A very interesting theoretical question is whether there is a truly deconflned quantum 
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critical point at the transition from superconductivity to ferromagnetism. Our results are 
consistent with this, and they imply that if there were an ordered phase, it would appear only 
at a tiny scale: in the RG flow, the growth of different terms in the interaction competes with 
the suppression by self-energy effects at the Van Hove points. At ti ~ 0.34, the growth 
tendencies of d-wave pairing and ferromagnetic correlations cancel one another, leaving 
the suppression of the quasi-particle weight as the dominant effect, which drives a drastic 
downturn of the stopping scale at that point and, in absence of further competition, will 
suppress that scale to zero. One possible shielding of the QCP would be the appearance of 
a c?-wave phase which is exclusively driven by the vicinity of the zone diagonals. At very 
weak coupling U , it could appear only below scales of the order exp(— ■^). We have also not 
seen it at the values of U and in the scale regime we study. A more detailed investigation 
using the effective action we derived here is under way. 
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